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FOREWORD 


This  report  describes  work  accomplished  during  1980  and  1981  in  the  program. 
Experimental  Mechanics.  This  program  is  funded  by  the  Lockheed  Independent 
Research  and  Development  Program,  and  administered  by  the  Structures  and 
Materials  Laboratory. 

Acknowledgement  is  due  K.  N.  Lauraitis  and  J.  T.  Ryder  for  their  helpful 
discussions,  and  to  R.  C.  Young,  F.  L.  Imera,  R.  LaForce  and  F.  M.  Pickel  for 
special  care  in  specimen  fabrication  and  testing. 
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ABSTRACT 


Buckling  and  crippling  of  a  structural  element  is  governed  by  its  elastic 
properties.  Nevertheless,  a  complete  definition  of  the  elastic  constants  of  a 
graphite/epoxy  laminate,  including  the  transverse  shear  moduli  (important 
under  high  compressive  stress),  has  never  been  accomplished.  For  this  pur¬ 
pose,  tension,  compression,  and  bending  tests  were  conducted  on  square-section 
specimens  cut  0-degrees  and  90-degrees  to  the  fiber  axis  from  64-ply  unidirec¬ 
tional  T300/5208  laminate.  Axial  test  specimens  were  instrumented  with  tee 
gages  to  obtain  data  on  extensional  moduli  and  Poisson's  ratios.  Bending 
deflections  were  measured  with  high  p-ecision  and  analyzed  using  a  rigorous 
finite  difference  solution  for  stress  distribution  to  derive  the  shear  moduli 
along  with  the  extensional  moduli  acting  in  bending.  Redundant  determinations 
were  obtained  for  the  nine  elastic  constants  of  Hooke's  law  as  generalized  for 
specially  orthotropic  material. 

Extensional  moduli  obtained  under  axial  load  and  in  bending  agreed  with  each 
other  to  within  one  percent.  A  self-consistent  set  of  properties  was  deter¬ 
mined  which  were  generally  in  accordance  with  published  data,  and  which 
established  the  material  as  being  transvers  .y  isotropic  and  conforming  to  the 
linear  theory  of  elasticity,.  The  value  obtained  for  the  shear  modulus  G 12 
(=G^)  was  4500  MPa  (0.65  Msi),  some  twenty  percent  lower  than  that  generally 
accepted  on  the  basis  of  445°  tension  test  data.  The  value  of  the  transverse 
shear  modulus  G ^  was  3650  MPa  (0.53  Msi),  in  accordance  with  published  re¬ 
sults  of  ultrasonic  tests. 
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INTRODUCTION 


In  the  design  and  analysis  of  metal  aircraft  structures  which  buckle  under 
compressive  loads,  and  those  which  deform  under  transverse  bending  load,  it  is 
seldom  necessary  to  consider  the  effects  of  shearing  deformation.  The  trans¬ 
verse  shear  stiffness  of  graphite-epoxy  composite  laminate,  however,  is  no 
more  than  a  fifth  that  of  aluminum  alloy,  and  the  adverse  effects  of  shearing 
deformation  will  require  consideration  in  many  practical  applications  of 
advanced  design.  Reference  1,  for  example,  indicates  a  reduction  in  com¬ 
pressive  buckling  stress  of  fifty  percent  for  stiffener  elements  designed  for 
high  stress  (width  to  thickness  ratio  of  10).  An  evaluation  of  the  transverse 
shear  modulus  of  this  material  to  a  reasonably  high  confidence  level  is  there¬ 
fore  necessary  for  applying  it  efficiently  in  aircraft  construction. 

The  elastic  properties  of  composite  laminate  are  usually  represented  as 
specially  orthotropic.  In  this  analytic  model,  unidirectional  laminate  has 
three  shear  moduli:  the  in-plane  shear  modulus  G12,  and  two  transverse  moduli 
and  ^23’  There  are  several  direct,  relatively  simple  test  methods  for 

evaluating  but  very  few  for  the  transverse  moduli.  A  selected  biblio¬ 

graphy  is  presented  in  References  2  through  15. 

The  most  commonly  used  method  for  evaluating  the  in-plane  shear  modulus  of 
composite  laminate  is  the  45°  tensile  test  (Reference  5).  In  this  test,  a 
tee  or  rosette  strain  gage  provides  bi-directional  stress-strain  data;  two- 
dimensional  transformation  of  stresses  and  strains  identifies  shear  stress  vs. 
shear  strain  in  the  plane  of  the  laminate.  While  this  test  is  concerned  only 
with  response  in  the  plane  of  the  laminate,  the  material  is  usually  assumed  to 
be  transversely  isotropic  on  the  basis  of  its  structure,  in  which  case  G^ 
=C12.  Investigations  have  demonstrated  this  method  to  be  equivalent  (probably 
to  the  limitations  of  experimental  control)  to  the  cross-sandwich  beam,  thin 
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tub®  in  torsion,  ten-deg  off  axis,  panel  shear,  and  rail  shear  tests  (Refer¬ 
ences  3,  6,  8).  Tests  of  T300/5208  graphite/epoxy  of  0.62-0.67  fiber  volume 
fraction  by  the  +95°  tension  method  have  indicated  =  0 .7**  -  0.80  Msi 

(References  9.18).  This  value  is  generally  understood  to  be  the  initial  slope 
of  the  stress-strain  relation.  A  possibly  significant  aspect  of  all  the  test 
methods  in  this  category  is  that  they  usually  furnish  decidedly  non-linear 
shear  stress-strain  curves  (for  example,  see  References  3,  6,  8,  9,  and  16). 

The  complete  set  of  material  elastic  constants  can  be  evaluated  from  the  speed 
of  acoustic  wave  propagation  in  ultrasonic  tests  of  specimens  of  special 

geometries.  Such  experiments  indicate  G12  and  to  be  about  thirty  percent 
higher  than  obtained  with  the  +  A5°  static  tension  test;  Reference  7,  for 

example,  gives  1.03  Msi  for  this  property,  with  G ^  -  0.527  Msi.  Engineers 

hesitate  to  use  such  values  in  design  because  a  one-to-one  relationship 

between  microscopic  dynamic  oscillations  and  macroscopic  static  deformation 
has  not  been  demonstrated  for  graphite-epoxy;  indeed,  even  in  isotropic 

metals,  the  two  techniques  usually  provide  slightly  different  results. 

Since  one  of  the  engineering  uses  of  the  shear  modulus  is  to  improve  the 
analytic  prediction  of  buckling  phenomena,  a  practical  approach  is  to  test 
simple  buckling  specimens,  determine  the  reductions  in  critical  load  due  to 
shear  effects,  and  back-figure  to  derive  applicable  values  of  the  shear 

moduli.  This  method  was  used  to  analyze  some  1800  column  test  data  in 

Reference  9.  For  T300/5208  laminate,  this  approach  indicated  room  temperature 
values  of  G^  =  0,30  Msi  and  G^  =  0.0M5  Msi.  Besides  being  unexpectedly  low, 
these  values  conflict  with  the  assumption  of  transverse  isotropy. 

In  the  past  few  years  the  Air  Force  Materials  Laboratory  has  developed  a 

direct  test  method  for  evaluating  two  of  the  shear  moduli.  This  method 

utilizes  a  half-inch  diameter,  six-inch  long  cyclinder  cut  90  degrees  to  the 

fiber  direction  from  thick  unidirectional  laminate.  Shear  strain  gages  are 

mounted  on  the  surface  at  0  degrees  and  90  degrees  to  the  fiber  and  the 

12. 

cylinder  is  loaded  in  torsion  (Reference  f>~) .  Tests  of  T300/5208  of  fiber 
volume  fraction  0.63  reported  in  Reference  j^,agive  values  of  0.866  and  0.506 
Msi  for  G.J2  and  G^,  respectively. 
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Deserving  of  mention  for  its  simplicity  and  low  cost,  a  method  known  as  the 
Iosipescu  shear  test  has  recently  been  applied  to  composite  laminates 
(Reference  15).  This  test  uses  a  specimen  only  slightly  larger  and  harder  to 
fabricate  than  that  for  the  ASTM  D  2344  short  beam  shear  test,  and  it  makes 
material  characterization  possible  for  both  shear  stiffness  and  strength  in 
all  three  planes.  The  method  must  be  confirmed  by  test  results  obtained  under 
loading  conditions  realistic  for  structural  applications,  since  (as  with  most 
methods)  the  geometry  is  unusual  and  the  purity  of  the  stress  condition  might 
be  questioned. 

A  method  that  is  appealing  from  a  practical  standpoint  determines  the  shear 
stiffness  from  simple  flexure  tests  in  accordance  with  accepted  handbook-type 
formulas  for  bending  and  shear  deflection  (References  10  and  14).  Application 
of  this  method  in  Reference  10  indicated  of  16-ply  T300/5208  laminate  to 
be  0.40  Msi.  This  method  suffers  from  the  crudeness  of  the  handbook  approxi¬ 
mation  made  for  shear  deflection,  as  well  as  from  the  difficulty  of  achieving 
adequate  experimental  precision. 

The  basis  for  an  accurate  analysis  of  flexure  test  data  was  established,  in  j 

Reference  19,  by  the  development  of  a  computer-aided  solution  for  the  plane 
stress  representation  of  a  simple  3-point  beam  of  orthotropie  material.  This 
solution  provides  a  precise  determination  of  the  complete  normal  and  shear 
stress  distribution,  and  so  permits  accurate  prediction  of  the  bending  and  the 
shearing  components  of  elastic  deflection.  This  predictive  calculation  can  be 
applied  to  determine  the  values  of  the  elastic  moduli  which  best  fit  redundant 
sets  of  deflection  data  obtained  in  flexure  tests  at  different  spans.  The  in¬ 
vestigation  described  in  this  report  was  planned  to  make  use  of  this  approach, 
in  order  to  evaluate  the  elastic  constants  of  graphite/epoxy  laminate  under 
loading  conditions  representative  of  typical  structural  applications. 
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OBJECTIVE 


The  express  purpose  of  this  investigation  was  to  evaluate  the  transverse 
shearing  stiffness  of  T300/5208  unidirectional  laminate  when  subjected  to 
simple  bending. 

In  the  completion  of  this  objective,  answers  to  the  following  additional 
questions  were  also  sought: 

o  Which  if  any  of  the  existing  test  methods  for  shear  modulus 
could  be  confirmed. 

o  What  causes  the  disagreement  among  values  obtained  by  different 
test  methods. 

o  Can  the  material  be  considered  transversely  isotropic  or  not. 

o  Is  more  fundamental  study  required  on  the  nature  of  the  material 
and  the  elastic  model  used  to  represent  its  response. 
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APPROACH 


Four  elastic  moduli  govern  the  plane  stress  or  plane  strain  response  of 
orthotropic  material.  In  the  x-z  plane  these  are  the  two  extensional  moduli 
E„„  and  E__,  the  transverse  shear  modulus  Gv_,  and  the  Poisson’s  ratio  v  . 

XX  zz  xz  xz 

The  bending  stiffness  of  a  flexure  specimen  is  a  strong  function  of  Exx,  a 
weak  function  of  Gv_,  and  almost  independent  of  E__  and  v  .  Therefore,  two 
flexure  tests  obtained  on  the  same  beam  at  different  spans  provide  sufficient 
information  to  evaluate  Exx  and  Gxz  to  good  precision.  To  do  this,  however, 
it  is  necessary  to  relate  observed  deflections  to  the  moduli  through  an 
accurate  representation  of  the  stress  distribution. 

The  two-dimensional  stress  distribution  for  an  orthotropic  beam  can  be  solved 

quickly  and  to  any  desired  degree  of  accuracy  by  the  computer-aided  relaxation 

19  , 

technique  described  in  Reference  p.  This  is  a  rigorous  finite  difference 
analysis  which  derives  values  of  the  stress  function  over  a  matrix  of  interior 
points  to  suit  a  specified  set  of  boundary  conditions.  Once  the  stress 
function  matrix  is  determined,  stresses  are  readily  computed,  and  the  strains 
and  deflections  follow  upon  selection  of  the  moduli. 

The  investigation  was  planned  to  take  advantage  of  conventional  axial  load 
tests  to  identify  the  extensional  moduli  and  the  Poisson’s  ratios  of  a  sample 
unidirectional  laminate.  Identical  square  section  specimens  cut  at  0  degrees 
and  90  degrees  to  the  fiber  axis  from  64-ply  graphite-epoxy  were  used  for  all 
tests.  A  number  were  tested  under  axial  load;  these  were  instrumented  with 
strain  gages.  Replicate  specimens  were  tested  in  three-point  bending,  the 
stiffness  being  determined  by  measurement  of  central  deflection.  Each  bending 
specimen  provided  data  in  two  planes  of  loading  under  identical  test 
conditions.  The  unreduced  test  data  therefore  furnished  qualitative 
information  regarding  transverse  isotropy  and  the  relative  magnitudes  of  G^ 
and  02y  The  experimental  approach  is  illustrated  schematically  in  Figure  1. 
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Deflection  measurements  made  In  the  flexure  tests  included  deformation 
occurring  in  the  test  fixture  and  Hertzian-like  bearing  deformation  of  the 
surface  of  the  beam  specimen.  These  effects  comprised  an  error  that  was 
significant  for  tests  at  very  short  span.  To  provide  correctional  data, 
supplementary  tests  were  conducted  on  specimens  loaded  as  a  block  in  bearing, 
and  on  the  test  fixture  with  the  specimen  removed.  The  correction  for  bearing 
deformation  was  derived  from  the  block-bearing  test  data  with  the  aid  of  two 
additional  computer-aided  relaxation  analyses,  by  which  the  Hertzian 
deformation  occurring  in  the  block  specimen  was  related  to  that  occurring  in 
the  beam. 

Reduction  of  the  flexure  test  data  was  achieved  by  applying  the  analysis  for 
beam  deflection  to  each  test  span  and  condition,  using  a  trial-and-correction 
procedure  to  obtain  values  of  the  moduli  which  best  fit  the  multiple  test 
points.  The  procedure  is  illustrated  by  the  the  flow  chart  of  Figure  2. 
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SPECIMENS 


To  furnish  stock  for  test  specimens,  a  64-ply  unidirectional  panel  was 
fabricated  of  T300/5208  low-resin  content  (35  t)  prepreg  tape.  The  layup  with 
symmetrical  bleeder  plies  was  placed  under  vacuum  in  a  press  for  the  standard 
cure  cycle  (45  minutes  at  132°C  followed  by  2  hours  at  179°C);  no  post-cure 
was  applied. 

To  remove  resin-rich  surface  material,  a  grinding  operation  using  water 
coolant  was  performed  on  both  surfaces  of  the  panel.  Square  specimens  of 
ten-inch  length  were  then  cut  0  degrees  and  90  degrees  to  the  laminate  axis. 
Final  grinding  reduced  all  sections  to  a  uniform  10.59  x  10.59  ran  (0.269  x 
0.269  inch).  Section  dimensions  were  measured  to  + 0.0001  inch  at  stations 
spaced  every  inch  along  the  length. 

Photomicrography  of  two  specimens  after  test  showed  a  high  degree  of 
homogeneity,  although  traces  of  the  ply  boundaries  indicated  irregular 
distortions  from  the  plane.  Photomicrographs  are  presented  in  the  Appendix. 

Resin  analysis  performed  per  ASTM  D3171  on  samples  taken  from  the  finished 
specimens  indicated  a  laminate  specific  gravity  of  1.61  and  a  fiber  volume 
fraction  of  68.6  percent.  The  specimens  were  dried  at  66  °C  O50°F)  in  vacuo 
for  two  months  prior  to  test,  a  procedure  estimated  to  reduce  the  moisture 
content  to  less  than  0.1  percent. 

Three  specimens  of  each  orientation  were  used  for  axial  tests  to  obtain 
extensional  properties  at  0  degrees  and  90  degrees  to  the  fiber  direction. 
The  problem  of  applying  test  load  in  the  thickness  (3)  direction  was  met  by 
cutting  short  lengths  from  specimens  tested  under  axial  load  and  reassembling 
the  short  sections  by  adhesive  bonding  to  form  3-tier  and  5-tier  compression 
specimens. 
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Additionally,  to  obtain  specimens  for  axial  load  in  the  3-direetion  which 
were  geometrically  similar  to  the  others,  cubical  elements  were  cut  from  a 
fourth  90-degree  specimen,  reoriented  through  90  degrees,  and  bond-sssembled 
as  indicated  in  Figure  1.  These  assembled  specimens  were  intended  for  test 
under  axial  tension,  but  bond  failure  prevented  tension  loading  of  useful 
magnitude.  One  of  these  specimens,  which  incorporated  seventeen  cubical 
elements,  was  salvaged  for  axial  compression  testing  by  squaring  the  ends. 

Three  additional  specimens  each,  of  0°  and  90°  orientation,  were  selected  for 
flexure  tests.  When  supported  for  beam  tests  the  central  load  was  always 
applied  at  the  specimen  midstation,  with  excess  specimen  length  extending 
symmetrically  beyond  the  beam  supports.  Maximum  variation  in  measured 
section  dimensions  over  the  test  span  for  any  of  these  specimens  was  less  than 
0.3  percent. 
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EXPERIMENTAL  PROCEDURE 


Axial  Tests 


Axial  test  specimens  were  instrumented  with  tee  foil  gages  centrally  mounted 
on  all  four  faces  and  tested  in  a  universal  testing  machine.  Gage  output  vs. 
load  data  were  automatically  recorded  and  reduced  with  the  Rye  Canyon 
centralized  digital  data  acquisition  system,  which  has  a  demonstrated  overall 
accuracy  better  than  0.05%  of  selected  range.  Tension  loading  not  exceeding 
0.003  strain  was  first  applied  using  Templin-type  mechanical  grips.  The 
specimens  were  then  reduced  to  approximately  10  cm  (9  inch)  length,  the  ends 
squared,  and  compression  loads  applied  between  flat,  adjustable  platens.  In 
compression  tests  of  built-up  specimens  in  the  thickness  (3)  direction,  the 
element  bearing  the  gages  was  always  the  central  member  of  a  stack  of 
identical  and  identically  oriented  elements,  in  order  to  minimize  end  effects. 

Representative  stress-strain  curves  obtained  in  these  tests  are  presented  in 
the  Appendix. 

Flexural  Tests 


Flexural  tests  were  conducted  in  the  special  fixture  shown  in  Figure  3-  toad 
was  applied  and  reacted  through  hard  steel  surfaces  of  3.175  nan  (0.125  inch) 
radius.  The  width  of  the  loading  nose  matched  that  of  the  specimen. 
Deflection  was  measured  between  the  loading  nose  and  the  fixture  base  with  a 
clip  gage  (MTS  632.02B01).  The  Rye  Canyon  centralized  digital  data  system  was 
used  to  acquire,  record,  and  reduce  load-deflec  ion  data.  End-to-end 
calibration  of  the  deflection  measurement  and  recording  system  against  a 
certified  Templin  extensometer  calibrator  indicated  overall  accuracy  and 
linearity  of  ♦  1.3  micrometer  (0.00005  inch)  maximian  deviation  over  a  range  of 
2.5  mm  (0.10  inch) . 
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Figure  3  -  Flexure  Test  Specimen  in  Test  Fixture 
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Tests  were  conducted  at  eight  different  spans  on  each  specimen,  providing 
replicate  data  from  three  specimens  for  each  of  the  four  beam  orientations. 
Test  span  lengths  had  to  be  short  enough  so  that  shear  deformation  was  large 
compared  with  the  expected  error  in  measurement,  yet  not  so  short  that  the 
error  caused  by  local  bearing  deformation  became  intolerable.  Maximum  applied 
loading  corresponded  to  a  beam  theory  strain  Mc/EI  of  about  0.004. 
Representative  deflection  vs.  load  curves  obtained  in  these  tests  are 
reproduced  in  the  Appendix.  Except  for  initial  softness,  which  may  be 
attributed  to  bearing  on  asperities  and  possibly  to  slight  accommodation  in 
torsion,  load-deflection  relations  were  quite  linear. 

The  measured  load  vs.  deflection  included  load  carried  by  the  clip  gage  as 
well  as  that  on  the  test  specimen.  A  test  was  conducted  on  the  clip  gage 
alone,  measuring  deflection  with  a  dial  gage,  to  evaluate  this  quantity.  The 
results,  which  are  included  in  the  Appendix,  provide  a  correction  which  was 
important  for  the  longest  test  spans. 

A  supplementary  test  was  performed  to  evaluate  the  error  introduced  into  the 
flexure  test  data  by  deformations  which  occurred  in  components  of  the  test 
fixture  included  between  the  gage  points.  For  this  test  the  loading  nose  was 
brought  into  contact  with  the  base  of  the  fixture,  and  the  load-deflection 
characteristics  recorded  as  before.  The  resulting  curve  is  reproduced  in  the 
Appendix.  (These  measurements  provided  only  an  approximation  to  the  behavior 
in  the  flexure  test  because  some  deformation  occurred  in  the  base  under  the 
loading  nose  instead  of  at  the  beam  supports). 

Tests  were  also  conducted  with  load  applied  to  representative  specimens  as  in 
the  flexure  tests  but  with  the  specimens  resting  directly  on  the  flat  base  of 
the  test  fixture,  in  order  measure  specimen  bearing  deformation.  Deflection 
vs.  load  curves  obtained  in  these  tests  are  included  in  the  Appendix. 
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TEST  RESULTS  AND  ANALYSIS 


Axial  Tests 


Compliance  data  under  axial  stress  were  derived  by  manually  evaluating  the 
slopes  of  the  stress-strain  curves  furnished  by  the  data  processing  system. 
This  is  an  operation  that  produces  repeatable  values  within  one  half  to  one 
percent  if  the  slopes  are  truly  linear.  Linearity  existed  over  the  major 
portion  of  the  curve  for  all  tests  except  those  for  tension  in  the  fiber 
direction,  in  which  stiffness  increased  slightly  with  load.  In  these  cases 
secant  data  were  obtained  over  several  ranges. 

Corrections  for  strain  gage  transverse  sensitivity  were  applied  after 
evaluation  of  the  slopes.  Tabulations  of  the  material  compliance 
coefficients,  obtained  by  averaging  data  from  back-to-back  strain  gages,  are 
recorded  in  the  Appendix.  Mean  values  and  coefficients  of  variation  are 
presented  in  Table  1. 

With  only  moderate  exception,  these  data  are  self-consistent  and  fit  the 
usually  assuned  model  which  displays  symmetric  coupling  effects  (S^  = 

S.J2  =  i ,  S23  =  S^2)  and  transverse  isotropy  ($22  =  S33’  S12  =  S13^*  *  set 
of  coefficients  for  this  model,  obtained  by  averaging  all  determinations 
bearing  on  any  particular  value,  is  as  follows  (units  of  ye/MPa  (ye/psi)): 


wiJ 

i.j  =  1.2.3 


6.92  (0.048) 
-2.2  (0.015) 
-2.2  (0.015) 


-  2.2  (0.015) 
94  (0.65) 

-  46  (0.32) 


-  2.2  (0.015) 

-  46  (0.32) 
94  (0.65) 


Corresponding  values  of  the  engineeering  constants  are  listed  in  Table  2. 
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TABLE  2  RESULTS  OF  AXIAL  LOAD  TESTS 


n 

Value 

Coef 

Var'n 

% 

Extensional  Moduli: 

E11 

14 

21.0  Ms! 

2.1 

E22'  E33 

27 

1.54  Msl 

2.9 

Poisson's  Ratios: 

V12’  v13 

14 

0.31 

3.8 

v23’  V32 

13 

0.49 

7.3 

V21*  v3t  (=  Vt2E22/E.1> 

13 

0.023 

. 

8.8 

'm&Lockheed 
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Flexure  Tests 

Recorded  test  results  were  provided  by  computer  in  the  form  of  load-deflection 
curves,  which  were  analyzed  manually  to  obtain  values  of  slope.  The  resultant 
deflection  rate  data  (yin. /lb)  and  the  corresponding  range  of  loading  over 
which  the  behavior  was  essentially  linear  are  listed  in  Table  3. 

The  correctional  data  also  were  reduced  to  terms  of  deflection  rate.  The 
correction  for  clip  gage  load  was  incorporated  by  combining  it  with  the 
average  deflection  rate  obtained  from  the  three  specimens  tested  under  each  of 
the  thirty-two  test  conditions.  These  results  are  included  in  Table  3. 
Bearing  deformation  effects  were  non-linear  in  nature,  and  reduction  of  the 
data  produced  the  deflection  rate  vs.  load  relations  presented  in  Figure  4.  A 
preliminary  analysis  of  the  flexure  test  data  was  required  before  the 
correction  for  bearing  deformation  could  be  applied. 

The  finite  difference  relaxation  technique  of  Reference  19,  which  was  used  in 
analyzing  the  flexure  test  data,  is  based  on  the  stress  function  relation  for 
orthotropic  material: 


TABLE  3.  RESULTS  OF  BENDING  TESTS 


DEFLECTION  RATE  (10'*  IN./LB.) 


LR  29763 


LOAD  P  (LB) 

Figure  4  -  Experimental  Deflection  Rates  (Slopes 
of  Recorded  Load  -  Deflection  Curves) 


Lockheed 

-California  Company 


LR  29763 


The  stress  function  solution  was  applied  to  treat  the  central  portion  of  the 
test  specimen  out  to  a  distance  2h  (twice  the  beam  depth)  from  the  point  of 
load  application.  Deflection  of  the  remainder  of  the  span  was  calculated  by 
beam  theory.  The  net  of  point  values  was  spaced  h/10  in  the  spanwise  (x) 
direction  and  h/20  in  the  depth  (z)  direction.  The  model  is  shown 
schematically  in  Figure  5.  Preliminary  studies  showed  that  this  analysis 
provided  the  desired  precision  in  deflection  calculation. 

Preliminary  solutions  of  the  stress  function  matrix  were  derived  for  each  of 
the  flexure  test  conditions,  using  elastic  constants  determined  in  the  axial 
tests  and  assuming  values  of  the  shear  moduli.  The  case  of  beam  load  in  the 
2-3  plane  (obtained  in  the  90-degree  beam  with  fibers  oriented  transverse  to 
the  load)  was  analyzed  as  plane  strain,  the  others  as  plane  stress.  Best-fit 
values  of  Exx  and  Gxz  were  then  derived  by  placing 

(dw  /dP)  =  a  (dw,  /dP).  +  b(dw  /dP),  (4) 
test  i  b  i  si 


Here  dw. /dP  and  dw  /dP  are  deflection  rates  attributable  to  bending  stress  and 
b  s 

to  shearing  stress,  respectively,  calculated  at  the  beam  midplane;  a  and  b 
are  error  factors.  Defining,  for  simplification. 


dw  /dP 
s 

dw,  /dP 

D 


dw  /dP 
test 

dw,  /dP 
b 


(5) 


Then 


yA  =  a  +  bxA  (6) 

and  a  and  b  may  be  derived  by  linear  regression  of  y  on  x.  If  bending  and 

shear  deformation  are  assumed  to  occur  independently  of  each  other,  a  and  b 

become  correction  factors  for  the  moduli  E  and  G  .  The  error  in  this 

XX  xz 

asssumption  becomes  negligible  if  the  correction  is  small,  making  iterative 
improvement  of  the  relaxation  solution  a  practical  measure. 
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Correction  of  the  flexure  test  data  was  then  made  for  bearing  deformations 
occurring  between  the  surface  at  which  load  was  applied  and  the  midplane  of 
the  beam.  Analysis  was  required  to  derive  this  correction  quantity  from 
measurements  taken  surface-to-surface  on  a  block  which  carried  no  bending 
stress.  General  solutions  for  the  two  cases  (viz.,  block  in  bearing  supported 
at  the  opposite  surface,  and  section  of  a  beam  supported  by  moment  and  shear 
at  the  ends)  were  obtained  by  extensions  of  the  stress  function  relaxation 
analysis  described  above  which  used  closer  grid  spacing.  For  each  of  the 
flexure  test  conditions,  a  set  of  specific  solutions  for  deformation  under  the 
two  support  conditions  was  derived  for  a  range  of  values  of  the  width  over 
which  the  central  load  was  distributed.  The  block-bearing  test  data  could 
then  be  used  to  evaluate  a  load  distribution  parameter  for  each  condition.  By 
assuming  this  parameter  applied  also  to  the  beam-supported  bearing  problem, 
the  applicable  solution  was  identified  and  the  bearing  deformation  occurring 
between  the  loading  nose  and  the  beam  midplane  was  evaluated. 

This  analytic  procedure  is  considered  only  an  approximation,  which  is 
especially  limited  for  short  spans  and  high  loads,  but  adequate  as  a 
correction  procedure.  To  facilitate  the  correction  process,  the  analytic 
solutions  are  arranged  in  nomographic  fashion  in  Figure  6.  Entry  to  the  lower 
set  of  curves  with  the  block-bearing  deflection  rate,  obtained  from  Figure  4 
at  the  median  of  the  load  range  of  interest,  lead  via  the  upper  curves  to  the 
3urface-to-midplane  deflection  rate  in  the  flexure  test.  Data  derived  by  this 
method  for  the  loading  station,  and  (under  half  the  beam  load)  for  the  support 
station,  were  applied  to  correct  the  flexure  data  in  Table  3.  Also  included 
in  Table  3  are  corrections  for  deflection  occurring  in  the  fixture,  taken  from 
Figure  4. 

With  the  flexure  test  data  refined,  additional  iterations  of  the  solution  for 
Exx  and  Gxz  were  performed  until  best-fit  values  to  the  eight  test  points  of 
each  set  were  obtained.  Variation  in  the  Poisson's  ratios  from  the  values 
determined  in  the  extensional  tests  was  not  investigated,  other  than  to 
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demonstrate  the  effect  to  be  too  small  to  utilize.  The  results,  which  provide 
at  least  one  and  in  some  cases  two  determinations  of  the  extensional  and  the 
shear  moduli,  are  presented  in  Table  4  and  summarized  in  Table  5.  Extensional 
and  shear  moduli  derived  from  beams  of  different  fiber  orientation  are  seen  to 
be  in  good  agreement,  and  agreement  of  the  two  transverse  shear  moduli  G12  and 
(from  0°  and  90°  beams,  respectively)  again  confirms  the  material  as  being 
transversely  isotropic.  Furthermore,  the  transverse  shear  modulus  is  *n 
almost  exact  agreement  with  the  value  expected  by  transverse  isotropy,  namely, 

G23  =  E22/2(U  V23}  °r  *52  Msi" 

The  fit  to  the  test  data  achieved  by  iterative  refinement  of  the  relaxation 
solution  is  indicated  in  Figure  7.  While  this  plot  demonstrates  that  the 
computation  provides  an  exceptionally  good  representation  of  total  beam 
deflection,  it  is  misleading  with  regard  to  the  precision  of  the  shear  modulus 
determination.  This  process  is  better  characterized  by  the  plot  of  final 
values  of  the  deflection  rate  ratios  presented  in  Figure  8. 
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TABLE  5.  SUMMARY  OF  FLEXURE  TEST  RESULTS 
(Extensional  and  Shear  Moduli  in  Msl) 


Test 

Condition 

0°,  3-dir 

0°,  2-dir 

90°,  3-dir 

90°,  1-dir 

N 

22 

22 

20 

21 

E11 

21.26 

20.96 

- 

- 

E22 

- 

- 

1.527 

1.552 

E33 

- 

- 

- 

- 

6,2 

- 

- 

- 

0.677 

G13 

0.652 

0.649 

- 

- 

G23 

- 

- 

0.536 

_ 

t 
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Figure  7.  Comparison  of  Computed  Solutions  for 
Beam  Deflection  with  Test  Data 
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I 
I 

1  Approximate  Method  for  Shear  Calculations 

|  Shear  deflection  occurring  in  beams  is  usually  estimated  by  assigning  a  value 

to  the  Timoshenko  shear  coefficient  K  in  the  relation, 

I 

!  “total  '  JJiT  d*  dx  +  dX  <7> 


I 

I 

I 

! 

I 

I 

I 

I 


Here  the  first  term  is  the  Bernoulli-Euler  bending  deflection  found  in 
handbook  formulas;  it  recognizes  no  restraint  of  section  warping.  The  second 
term  is  a  "shear  correction"  which  includes  all  other  effects,  presumed  to 
cause  a  shearing  type  of  deflection.  K  is  generally  thought  to  be  a  section 
property  that  varies  only  with  section  shape;  various  sources  cite  values 
ranging  from  2/3  to  1.0.  Consideration  of  the  shear  stress  destribution  in 
beams  makes  it  clear  that  K  may  vary  over  the  span  (Reference  19). 
Therefore,  any  value  which  is  to  be  used  outside  of  the  integral  as  in  (7) 
must  be  specific  not  only  for  load  and  section  but  for  span. 

The  results  of  the  current  study  can  be  utilized  to  develop  rational  values  of 
K.  Such  data,  obtained  by  introducing  the  appropriate  values  from  Table  4 
into  (7)  and  solving  for  K,  are  presented  in  Figure  9.  (For  the  90°  beams 
loaded  in  the  3-direction,  the  value  of  E  was  taken  as  E22/(1-  V  2^  V.^)  in 
accordance  with  representation  as  plane  strain.) 

Reissner  (Reference  20),  and  Nair  and  Reissner  (Reference  21),  working  from 
fundamental  energy  considerations,  have  established  upper  and  lower  bounds  for 
the  shear  deflection  that  occurs  in  a  beam  which  does  not  deform  locally  at 
the  load  and  support  stations.  In  the  90-degree  beam  loaded  in  the  1- 
direction,  the  high  fiber  stiffness  produces  boundary  conditions  which  are 
closely  represented  by  Reissner’ s  model,  and  the  resulting  values  of  K,  as 
plotted  in  Figure  9.  are  seen  to  lie  close  to  the  prescribed  bounds  calculated 
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from  Reference  21.  When  the  load  is  transverse  to  the  fiber,  however,  local 
deformation  at  the  loading  station  is  significantly  greater,  and  consequently 
K  for  very  short  spans  is  increased  and  may  exceed  1.0. 

These  differences  in  the  value  of  K  are  produced  by  variation  in  the  way  shear 
is  introduced  into  the  beam,  interacting  with  section  constraint.  Shear 
stress  distribution  adjacent  to  the  loading  station  for  the  90-degree  beam,  L 
=  1.02  inch,  loaded  in  the  3-  direction,  is  contrasted  with  that  for  load  in 
the  1-  direction  in  Figure  10  (data  from  the  relaxation  solutions).  The  high 
peak  stress  and  distorted  distribution  associated  with  loading  transverse  to 
the  fibers  tends  to  produce  large  section  warping.  The  central  section  cannot 
warp,  however,  and  adjustment  to  this  condition  reduces  overall  beam  de¬ 
flection.  Thus,  short  beams  with  low  E  zz  may  appear  to  be  stiffer  in  "shear" 
than  those  for  which  E  zz  is  large. 

Accuracy  in  the  calculation  of  K  is  adversely  affected  by  limitations  in  ac¬ 
curacy  of  the  current  set  of  relaxation  solutions.  While  these  data  have  a 
precision  of  better  than  0.1  percent,  they  suffer  a  small  error  as  a  result  of 
the  finite  dimensions  of  the  grid.  Such  errors  have  a  greater  effect  on  K 
than  on  the  shear  modulus,  because  after  Gxz  has  been  evaluated,  the  relax¬ 
ation  solution  values  enter  once  again  in  the  term  of  Eq.  (7),  and  in 
this  case  without  the  benefit  of  multiple  points.  For  cases  important  in 
design,  further  investigation  is  desirable  to  confirm  the  results  and  to 
establish  a  broader  basis  for  extrapolation  to  other  materials. 
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Overall 


Where  data  obtained  in  the  axial  load  testa  and  the  flexure  tests  overlap, 
only  small  differences  are  seen  to  exist,  and  these  discrepancies  probably 
represent  the  accuracy  of  the  test  methods.  By  taking  rounded  values  which 
are  consistent  with  each  other  and  the  indications  of  transverse  isotropy,  a 

indepen<ten"Hy  ev»iu*tle<i 

complete  set  of  the - nine-  elastic  constants^is  obtained  and  presented  in  Table 

6. 


These  data  represent  the  overall  results  of  replicate  tests  and  various  test 
conditions.  Scatter  of  test  results  with  different  test  procedures  is  no  more 
than  a  few  percent. 

These  results  were  derived  by  assuming  linear,  orthotropic  elasticity,  an 
assumption  justified  by  the  structure  of  the  material  and  the  linearity 
exhibited  in  the  basic  load-deflection  data.  The  fact  that  the  final  results 
provide  a  complete,  self-consistent  set  of  properties  confirms  this  model  and 
appears  to  refute  suggestions  of  non-linear  elasticity  relations  found  in  the 
literature  (Reference  16  and  17.  for  example).  However,  the  scope  of  the 
current  program  with  regard  to  stress  range  and  condition  is  insufficient  to 
resolve  this  question.  For  the  type  of  loading  and  the  material  tested, 
simple  c'thotropy  provides  a  satisfactory  model. 
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TABLE  6  OVERALL  SUMMARY 


Elastic  Constants 

of  T300/5208  Laminate 

As  Tested 

PROPERTY 

VALUE 

(GPa) 

(Msi) 

I 

E1  1 

145. 

! 

21.0  ! 

E22  •  E33 

10.7 

1.55  j 

G12  ’  G1 3 

4.50 

i 

0.65 

G23 

3.70 

0.52 

v12  ,V1 3 

0.31 

v23  ,v32 

0.49 

V21  ,V31 

0.023 
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CONCLUSIONS 


A  complete  set  of  the  nine  elastic  constants  experimentally  determined  for  a  64- 
ply  T300/5208  unidirectional  laminate  is  presented  in  Table  6.  These  data  are 
based  on  multiple  independent  experimental  determinations,  with  the  overall 
means  of  the  test  values  modified  slightly  to  obtain  a  consistent  set. 

The  results  agree  with  values  cited  in  the  literature,  except  in  the  case  of  the 
in-plane  shear  modulus  G12.  The  value  obtained  here  (0.65Msi)  is  somewhat  lower 
than  that  expected  from  +  45°  tension  tests  (0.80  Msi)  and  the  Pagano  torsion 
test  (0.87  Msi),  and  substantially  lower  than  reported  for  ultrasonic  tests 
(1.03  Msi) . 

The  results  obtained  confirm,  to  the  limits  of  the  experiment,  that  the  elastic 
stress-strain  relations  are  linear  and  that  the  material  is  transversely  iso¬ 
tropic. 

Calculations  of  the  Timoshenko  shear  coefficient  indicate  a  substantial  varia¬ 
tion  with  the  orientation  of  the  orthotropic  axes  of  the  beam  and  with  the  span. 
In  some  cases  the  calculated  values  exceed  by  a  considerable  amount  the  recog¬ 
nized  bounds  for  a  cantilever  which  has  slightly  different  end  conditions. 
Boundary  conditions  are  thus  seen  to  have  an  unexpectedly  strong  influence  on 
the  deformation  attributable  to  shear.  Such  an  effect  could  explain  variations 
in  performance  associated  with  different  test  procedures,  and  would  increase  the 
problems  of  analysis  and  design. 
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Figure  A- 1  -  Sections  of  Flexure  Test  Specimens  0°  No.  1  (left;  am'  yQ°  No.  1  (right),  at 
Magnification.  Damage  to  external  surface  of  0°  specimen  was  caused  by  test 
machine  grips. 
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table-  a-/.  compliance  coefficients 


TENSION  TESTS  IN  FIBER. 


10-40 


40-  SO 


COMPLIANCE  3 =  b  C;  J&Cj  Qjt/pst) 


SPECIMEN 


O-  tO  O  No.  t 


O  No.  2 


0°  No.  3 


O  No.  I  // 672/2.  H 3  0.04  834 

2*4  0.047  3  5 

O*  No.  z  mcl/t  NS  0.048  1  7 

2*4  0.04638 


0°  No.  3  /204t>/4  1*3  O  .047  H 

2*4  0  .0  4674 


O a  /Vo,/  N672j2  l* 3  0.04  65  7 

2*4  0  .  04  553 

0°  No. 2  11767/1  1*3.  0.  04648 

2*4  0.  04584 


10-40  A/ 


ME- AN 

COEF .  VAP  N 


0.0  46! 


- 0.0/50 
0.027 


^jLockheed 

-  California  Company 


-O.  0/525 


-0.014  7 
0.030 


LR  29763 


TABLE-  A-2  .  COMPLIANCE  COEFFICIENTS 

COMPR.ES. S ion [  T£S rs  /A/  r/BEJL  X>/RECT/ON  (t-  D/P.) 


SPECIMEN 

TEST 

STRESS 

sase 

COMPLIANCE  sr=  <5e;/z<r  (u€  /psl 

v  j  y  j 

MO. 

PAUSE 

(ksi) 

**3/ 

0°  No.  2 

7/8S//4 

B 

!+3 

2*4 

<J  .  <3435-2, 

0. 04846 

-0.01454 

-  0.0/46  5 

0°  Ho.  Z 

uasi/s 

/-/o 

M3 

0.  050/8 

-  0.01560 

2-*  4 

0.0430/ 

-  0.0/583 

0°  No.  3 

12046/2 

5--// 

/  +  3 

0.04351 

-  0 .0/530 

2^4 

0.048  78 

-0.0/330 

A/ 

6 

3 

3 

MEAN 

0.0432 

- 0.0/5 / 

-O.  O  146 

COEF. 

144 

0.  016 

0.036 

0 . 066 

'*&LocMteed 
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T A &L  E  A- 3.  COMPLIANCE  CQEFP/C/EUTS 

TENSION  TESTS  IN-PL AHE.  TRANSVERSE  f 2-D/R. 


SPECIMEN 


STRESS 

RAN&E 

(ksi) 


COMPLIANCE  S;:  -  be i  /b<Tj  juE/psi) 


30°  No.  /  7/67^2 


/ +  3  0.  6  486 
2*4  0,  6464 


-0.  6  1404 


- 0.3104 


30°  Nt i.  2  insi/l 
II76//2 
l/703/i 
11703/2 


O.  033 7 
O.  6  3/2 

O.  6S37 
Q.  6  575 

0.  0  3 60 
0.  0420 

0,  *357 

0.0  30  0 


-0.0/563 


- 0 .01553 


- 0.0/54  / 


-0.0/533 


-0.3050 


-0.3076 


-O  .3 043 


-0.3001 


30  No.  3  12053/1 


0.6385 
0.  641  7 


-0.0/527 


■0.3050 


-  0.0/56 
0.036 


■0.30  7 
0.007 


y^rLockheed 

CaHtom*  Company 
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TABLE  a- 6. 


COMPLIANCE  COEFFICIENTS 

COMPRESSION  TESTS  / N-PLANE ,  TRANSVERSE  ( 2-JD/ft 


SPECIMEN 


STKEii  COMPLIANCE  S;j  *  ie:  /btj 

PANte  *ASE 

(ksi) 


9 0° No.  I  /1676/S 

UC76,/<0 


0.  6276 
3.634/ 

0-6  4/3 
0.  6369 


■0.014  73 
-0.01373 


-0.30ZZ 

-0.2316 


SO  No.  2  /  2337/2 


0.6  IS* 
2*4-  O.  6/24 


-0.0!  S44 


-0.2967 


N 

(0 

3 

3 

MEAN 

0.628 

-0.  0/46 

-0.237 

COEF.  VAR.N 

0.01 9 

o.oss 

0.016 

'*&Lockheed 
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! 

I 


JA8LE  A-S .  COMPLIANCE  COEFFICIENTS 

I  COMPRESSION  TESTS  IN  THICKNESS  DIRECTION  (3-01*.} 


SPECIMEN 

\ 

TEST 

NO. 

STRESS 

RAN6£ 

(hsi) 

GfiCE 

compliance  s£j  -  ZeL/dcr-  (u6/p*i) 

533 

5/3 

SZ3 

0°  No.  2 

I1.032J5 

5  -10 

2+4 

0.4400 

-o.o  n  s 

(  5-TtEf£) 

I203zj(, 

4,-/3 

2+4 

0.4722 

-1 0.0/18 

30°  No.  2 

I20SS/2 

5-  9 

2  +  4 

0.  4433 

(3-TIEK& 

I120SE/3 

2-/2 

2  +  4 

0.4  4  58 

-0.3E9/ 

/20S8/5 

8-/6 

2+4 

0.4  4  02 

-0.3509 

30-90- Z 

I3054-) 4- 

1-3 

/  ■#  3 

0.4705 

-0.0/4  / 

(n-TiE$ 

2  +  4 

0.4878 

-0.3330 

I30S4/S 

2-4 

7  +  3 

0.4752 

-0.0/ 53 

2  +  4 

0.4704 

-0.3274 

N 

9 

4 

5* 

MEAN 

0.47/ 

-0.0/87 

-  0.  344 

COE  F.  VAR,'" 

0.0/4- 

0.072 

0.043 

'*&Lockheed 

-CaBfbmta  Company 
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